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Abstract 

The main objective of this paper is to study the efficiency of genus statistics in differ- 
entiating between different models of structure formation. Some of the models studied 
. ^ are over-simplified but aim at reproducing some of the features of the structure seeded 

by topological defects. We consider accretion onto static point masses that could ap- 
proximate accretion onto slow-moving cosmic string loops or other primordial point-like 
sources. Filamentary structures and wakes are also considered as an approximation to 

also made with predictions of genus statistics for Gaussian fiuctuations and with genus 
curves obtained from the CfA survey. A generic class of density models with wakes and 
filaments was seen to provide results comparable or better than Gaussian models for 
this suite of tests. 
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1 Introduction 

There are several models that aim at explaining the large scale structure we observe in 
the universe today. The two main paradigms for the origin of the seed perturbations are, 
on one hand, the inflationary scenario that produces Gaussian fluctuations and on the 
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other hand, topological defects produced at phase transitions in the early universe that 
produce a non-Gaussian spectrum of density perturbations, specially on small scales. 
While, in general, inflationary scenarios require great flne tuning of parameters in or- 
der for the required perturbation amplitude to be produced by the present time, models 
where the structure is seeded by topological defects depend mainly on a single parameter 
which usually is the energy scale of the phase transition. It turns out that for cosmic 
strings, global monopoles and textures the energy scale is that of GUT phase transi- 
tions. In a way these models seem to explain the origin of the seed perturbations more 
naturally than inflationary models and so they deserve close attention. At present most 
work is done on inflationary models and, because they produce Gaussian fluctuations, 
all statistics are completely described by the power spectrum; other statistics that could 
help to distinguish between inflationary and non-inflationary models are often neglected. 
Because both models produce a nearly scale-invariant spectrum of fluctuations and be- 
cause the uncertainties in the measurements are still large, the power spectrum is not 
the best statistic to distinguish between them. So, a higher order statistic is required. 
In this article we investigate the efliciency of genus statistics in distinguishing between 
different models of structure formation such as inflationary and cosmic string models 
and compare with observational genus curves obtained from the CfA survey. Some work 
on genus curves of isodensity contours for toy models of structure formation seeded by 
topological defects was previously done by Brandenberger, Kaplan & Ramsey (1993). 
By requiring that all the structures have the same size and mass they appear not to have 
properly taken into account the scaling solution. The size of the defect seeds increases 
with time (proportionally to the horizon) and the amplitude of the density perturbations 
induced by larger defects is smaller because they have less time to grow by gravitational 
instability. This effect is properly taken into account in the present paper and so we 
consider our work an improvement over that of Brandenberger et al. (1993). Robinson 
& Albrecht (1995) also performed a similar study with cosmic string wakes. Their cosmic 
string toy model consisted of a realization of a string power spectrum, where the phases 
of the Fourier modes were choosen at random, plus a single cosmic string wake. They 
concluded that the genus statistic is not a good discriminator between their model and a 
model without the wake included. As we will show in this paper the genus statistic is a 
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good discriminator between our cosmic string toy models and random phase models with 
the same power spectrum at least in the linear regime. So, we conclude that their model 
underestimates the presence of sheet-like features in the density field and, therefore, in 
this regard it should not be considered a good cosmic string toy model. 

This article is organized as follows. We start in section 2 by introducing the 
Zel'dovich aproximation, solving it for accretion onto static point masses, filaments and 
wakes assuming the dark matter to be cold. In the end of the section a modification 
to the Zel'dovich aproximation that accounts for the neutrino free streaming length is 
introduced when considering a hot dark matter model. In section 3 the genus statistic 
is described. Analytic results for Gaussian perturbations are given and topological mea- 
sures of departures from Gaussianity are introduced. In section 4 we describe the way 
we produce the fluctuations for the models considered and relate that to the results of 
section 2. In section 5 the results are presented. The dependence of the genus curves on 
the type and number of topological defects present is demonstrated and error bars due 
to sample variance are introduced. The density probability distribution for several of 
the models studied in this article are given and the parametrization of the genus curve 
discussed. The flnal section is a discussion of the results. 

In this article we shall assume that = 1 and h = iJo/(100 Km/sec/Mpc). 

2 The Zel'dovich approximation 
2.1 Cold dark matter 

We employ the Zel'dovich approximation (Zel'dovich 1970) to examine each model we 
use. We will assume the universe to be flat, with no cosmological constant and the dark 
matter to be cold. In the next section we will see how the Zel'dovich approximation can 
be modifled in order to describe hot dark matter as well. The equations for the evolution 
of the scale factor are the Friedman and Raychaudhuri equation that in the matter are 
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where py, is the background density. The position of the CDM particles is given by 

r = a(t)[q + ^(q,t)], (3) 

where q is the unperturbed comoving position of the particle and ?/> is the comoving 
displacement vector of the particle. In the presence of a perturbing string seed, the cold 
dark matter particle will obey 

dt 



2 Fgeed ~l~ Fjj^atter 7 (4) 



where the force Flatter due to the surrounding matter is given by 

F„,atter = -V$, (5) 

with the gravitational potential $ satisfying the Poisson equation 

= 47rGp„,atter = 47rGpb(l + S). (6) 

The fact that mass is conserved is described by the continuity equation that can be 
written in comoving coordinates to first order in ?/> as 

^+V,-(pa^$) = 0, (7) 

which can be integrated in the limit of small perturbations to give 

<^= — ~-Vq-^, (8) 
P 

again to first order in The Poisson equation can then be integrated to give 

Vr$ = ^^(r-3a^ll), (9) 
where ^/jy and ip± are defined by 

^ = ^||+^±, (10) 

VrX^||=0, (11) 



Vr-V'± = 0. (12) 

In the linear regime with (Itpl << |q|) we can use (4), (5) and (9), to write the Zel'dovich 
approximation. 

°(^+2^| + 35)j„=F.,,„, (13) 

/ 5^ a d \ ^ 

If the perturbations are irrotational to begin with, and if the source term is irrotational 
j_ = so that ?/> = ^/jij. Then we can write the Zel'dovich approximation as 

|-+2-|^ + 3-^^ = -Fseed =Sseed. (15) 

Ot"^ a ot a J a 
If the seed is a point mass Sgeed is given by 

__GM(q-cu) 

''pmass — i3 9 ' v / 

|q - qs| a-^ 

where qg is the comoving position of the point mass. If the seed is a line of mass then 
Sseed is given by 



__2GM4q-cu) 

'-'Imass — |2 o ' v / 

|q - qs| 

where qg is the comoving position of the point in the line of mass nearer to q and Ml 
is the mass per unit length of the line of mass. The Zel'dovich approximation can be 
solved using the Green's function method. For the case we are considering the Green's 
function is 



G(t,to) = t<to. (19) 

Using it we find the solution of the Zel'dovich equation for the static point mass and line 
of mass 



^»,— (x,<) = -^-5^ (-1 ln( -!--(- I . (2i: 



The formation of wakes can be modeled given an initial velocity of the form 




u^e(q), 



(22) 



where 




1, 



X • q > 0, 



(23) 



e(q) = 1, 



X • q < 0. 



(24) 



This is given to the particles in cosmic string models because of the conical spacetime 
with Ui = SttGiiVs^ where Vg is the string velocity. The solution of the Zel'dovich equation 
in this case is 



2.1 Hot dark matter 

If the dark matter is hot, then small-scale perturbations with wavelength smaller than 
the neutrino (or other hot dark matter candidate) free streaming length \ps will be 
erased. In order to properly describe the formation of structure in the context of a 
hot dark matter model one would need to solve the Boltzmann equation for the dark 
matter particles. However, it was shown by Perivolaropolous, Brandenberger & Stebbins 
(1990) in studying the clustering of neutrinos in cosmic string induced wakes, that most 
of the results can be described correctly using a naive modification of the Zel'dovich 
approximation. This modification is based on the fact that on average the dark matter 
particles will only start to collapse when the comoving free-streaming length has fallen 
below |q|. So we modify the Zel'dovich approximation in the context of a hot dark 
matter model by setting I?/;! = on scales |q| < \ps and evolving scales |q| > \ps as 
for the cold dark matter case. We consider a hot dark matter model with two neutrino 
species with sufficient mass to make Oj^ = 1. In particular we consider these two species 
to have the same mass nii^ = 46 eV. 

The comoving distance hot particles can move since t > teq is 




(25) 



\fs = 3v{t)tz{t), 
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where v[t) is the thermal velocity of the hot particles and z[t) is the redshift. The 
free streaming length grows in the radiation era attaining a maximum at teq and then 
decreases proportionally to t~^/^ in the matter era. The maximal free streaming length 
is 

^FS = SUeq^eq^^eq, (27) 

where 

v,^ = Tl'^lm, ^ 0.09, (28) 
so that the maximum free streaming length is given by 3/z~^Mpc. 

3 Genus statistics 
3.1 Definition of genus 

To measure the topology of isodensity contours we will use the Gauss-Bonnet theorem 
which relates the integrated Gaussian curvature (a local property) of a surface with the 
genus (a global property) of that surface. The Gaussian curvature of a two-dimensional 
surface at a particular point is 

K = (29) 
aia2 

where ai and a2 are the two principal radii of curvature at that point. A surface has a 
positive or negative Gaussian curvature respectively if the two radii of curvature point in 
the same or in opposite directions. For example, a sphere has a positive radii of curvature 
given by K = where r is the radius of the sphere. A cylinder has K = because 

one of the radii of curvature is infinite. Saddle points have negative curvature because 
ai and a2 have opposite signs. The Gauss-Bonnet theorem relates the integral of the 
Gaussian curvature over the surface with the genus in the following way 

J KdA = 47r(l -^), (30) 

where g is the genus of the surface and dA a surface element. The genus measures the 
number of closed curves that may be drawn on a surface without separating it. It can 
also be defined as 

g = number of compact regions — number of holes -|- 1. (31) 
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A curved surface may be approximated by a network of polygonal faces. When we use 
such a network to compute the genus we find that 



where Di = 2tt — ^ Vi is the angle deficit at a vertex in radians. In this case the curvature 
is effectively compressed into delta functions at the vertex. We used (32) to construct 
a numerical alghoritm in C in order to determine the genus of an isodensity surface 
applying the method sugested by Gott, Mellot & Dickinson (1986). The isodensity 
surface is constructed by binning the density onto a cubic lattice and identifying pixels 
with density above and below a certain threshold (Sc). Our program to compute the 
genus was tested against analytic results for Gaussian perturbations and also against 
well known topological configurations for which we knew the genus beforehand (e.g. a 
network of isolated cubes). 

3.2 Window function 

The data is smoothed with a Gaussian window function 

1 . . 



where Ag is the smoothing scale. This definition implies that the smoothing length is 



are always greater than the average interparticle spacing but not too large as to erase 
all the relevant features in the density map. The smoothing scales considered in this 
article were 6, 8, 10, 12, 16 and 20/z~"^Mpc in order to make a direct comparison with 
the results of Vogeley et al. (1994) for the CfA survey. 

3.3 Genus for Gaussian random fields 

A comparison with inflationary models predictions for the genus curve is essential if one 
has to decide which model describes better the kind of large scale structures we observe 
in the universe today. The genus curves of random flelds are well studied and some 
analytic formulae have been derived (Bardeen et al. 1986; Hamilton, Gott & Weinberg 
1986). The genus per unit volume is given by 







greater by a factor of \/2 than the usual width of a Gaussian. The smoothing scales 
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where v is the number of standard deviations above or below the mean density contour. 
The amphtude iV depends on the power spectrum of density fluctuations P(A;) as 

where 

If we smooth the structure on a scale Ag the power spectrum becomes 

P'(A;) = P(A;)e-'^'^^'/^ (37) 
and if the power spectrum is of the form P(A;) oc A;" then iV is given by 

Q I \ 3/2 



^ = ^ • 38 

(27r)'Aen 3 ; 

To compute the genus curve we must multiply (34) by the volume of the grid. In flg. 1 we 
plotted a random phase genus curve obtained for a P(A;) oc k power spectrum obtained 
respectivelly analytically, using (38) and (34), and numerically using our program to 
compute the genus. The two curves are almost identical as we should expect if the 
program to compute the genus is working properly. Small differences between the two 
curves can be atributed mainly to the choice of periodic boundary conditions. Smaller 
differences are also due to the sample variance (the volume of the box is not inflnite) 
and to small numerical imprecisions. The flgure is symmetric with respect to the vertical 
axis which puts in evidence the topological equivalence of positive and negative linear 
density perturbations for random phase models of structure formation. The shape of the 
random phase genus curve is independent of the power spectrum. For < 1 the genus 
is always positive, the surface has more holes than compact regions and so the surface is 
"spongelike". For > 1 the genus is negative and the surface has a lot of independent 
compact regions. For non-random phase distributions the genus curve will be, in general, 
asymmetric because the topological symmetry between high and low density regions is 
not expected in most cases. 
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3.4 Genus metastatistics 

To quantify departures from the random phase curve, Vogeley et ah (1994) used genus 
related statistics like the amplitude, the width, and the shift of the genus curve. The 
amplitude was defined as the amplitude of the best fit random phase curve. That was 
done by determining the amplitude of the genus curve that minimizes ■ The reason 
for the definition was that they wanted to compare the observational data from the CfA 
survey with predictions from Gaussian models; So, it seemed appropriate. Although, we 
want to make a comparison of present observations with non- Gaussian models we will 
retain the same definition in order to directly compare our results with those of Vogeley 
et al. (1994). 

The width of the genus peak W,j was defined as the difference between the zero 
crossings of the genus curve which for random phases is equal to W,j = — V- =2 
because, as we have seen before the genus for random phases is positive over the range 

— 1 < < 1 and negative elsewhere. This change of the genus sign is believed to 
coincide with the percolation thresholds for random-phase perturbations. In the range 

— 1 < < 1 both high ((^ > he) and low i^h < 5c) density phases percolate, while for 
u > 1 only the low density phase percolates and for u < —1 only the high density phase 
percolates. 

The last statistic they used was the shift Au of the genus curve which was quantified 
in the following form 

Au=^^ — , (39) 

/_i G{u)fitdu 

where G{u)^^^ is the measured genus curve and G{u)^^ is the best fit random phase cuve. 
A positive value of Au indicates a density distribution that is more 'meatball-like' (iso- 
lated cluster models) than random phase, while a negative value of Au is characteristic 
of a 'bubble-like' topology ('swiss chesse' topology). 

3.5 Parametrization of the genus curve 

For a Gaussian density field the volume fraction in the high density region is given by 

f{u >uo) = ^ I e-'"l^dt. (40) 



Vogeley et al. (1994) did not compute the mean and standard deviation u of the density 
distribution and express the genus curve as a function of u. Instead they determined 
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the genus curve as a function of the volume fraction in the high density region and 
used (40) to parametrize / as a function of v. Although for Gaussian perturbations 
these two ways of calculating the genus curves would give similar results, for some of 
the models we investigate in this article the genus curves are quite different due to their 
non-Gaussianity. 

3.6 Smoothing of the genus curve 

The genus curves were smoothed using a very simple procedure known as three-point 
boxcar smoothing (see for example Vogeley et al. 1994) that was shown to give better 
estimates of the true genus curve for Gaussian random phase models. It consists in 
determining the genus as 



where yi-\-i =1^^ + 0.1. 

4 Models and observations 
4.1 Toy models 

The toy models we study in this paper are simplified cosmic string models. It is a 
well-known fact that the gravitational effect of a slow-moving small loop can be well 
approximated by that of a static point mass, thus generating spherical accretion (see for 
example Vilenkin & Shellard 1994). Also a slow-moving wiggly string can be approx- 
imated by a line of mass, thus generating filamentary structures, while a fast moving 
string generates a wake. It is also known that cosmic string networks after the friction 
dominated era rapidly attain a scaling solution where the average properties of the net- 
work (such as the average number of defects and the average correlation length) remain 
the same at all times when scaled to the horizon size (Bennett & Bouchet 1990; Allen 
& Shellard 1990, Albrecht & Turok 1989). This means that although fluctuations laid 
down at later times are smaller in amplitude because they have less time to grow by 
gravitational instability, they will have a larger wavelength (in proportion to the horizon 
size). These facts are essential ingredients of the toy models we consider in this paper. 




(41) 
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Here we consider several toy models of cosmic string structure formation without 
taking into consideration the detailed properties of cosmic string networks known from 
numerical simulations. These models are respectively a network of fast-moving strings, 
a network of slow-moving wiggly strings and a network of slow-moving small loops. 
Both the filament and wake models should give a rough approximation to the structures 
seeded by cosmic strings. Although the network of slow moving small loops cannot be 
considered a realistic cosmic string toy model because loops of cosmic string produced 
by a cosmic string network are born with relativistic velocities and are in much higher 
number than what is considered in this article, there are scenarios which this toy model 
does approximate, notably those with loop nucleation during inflation. The kind of 
shapes we investigate also appear in other defect-seeded structure formation models like 
those seeded by global monopoles or global textures. Although we want primarily to 
test if the genus statistic is a good discriminator between different models of structure 
formation (specially between different non-Gaussian models) we also want to see if some 
of the features of these simplifled models match current observations. To properly test 
cosmic string models of structure formation one would need to go beyond this simplifled 
models and perform large-scale network simulations (Avelino & Shellard 1995). 

In flg 2 we plotted the isodensity contours for some of the CDM toy models con- 
sidered. We can see mainly fllamentary, wake-like, and spherical structures respectively 
in the 5 fllament, 5 wake and 25 sphere models. Other kinds of shapes can be obtained 
due to superposition of density perturbations generated by several defects. In the sphere 
model there are more small spheres than big ones because denser objects are generated 
later in this model when there are less defects inside the box. In the fllament and wake 
models we can see more smaller wakes and fllaments than larger ones. This is due to the 
fact that smaller objects are seeded earlier when there are more defects inside the box. 
These give rise to larger density perturbations and so to thicker objects in the density 
contour plot. 

4.2 The CfA survey 

Vogeley et al. (1994) studied the topology of large scale structure in the CfA Redshift 
Survey. This survey includes 12000 galaxies with limit magnitude mi, < 15.5 and 
allowed for the computation of the topology on smoothing scales from 6 to 20/z~"^Mpc. 
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They used genus statistics to test several variants of the cold dark matter cosmogony 
(CDM). All of them failed to match the observations to a high confidence level (> 97%) 
(Vogeley et al. 1994), even when evolution of the perturbations into the non-linear regime 
through the use of N-body codes was taken into account. This is a good motivation for 
our work based on non-Gaussian models of structure formation. Our results for the 
cosmic string toy models of structure formation are such that direct comparisons with 
the CfA genus curves are possible. 

4.3 Generation of fluctuations 

We apply the genus statistics to smoothing scales between 6/z~"^Mpc and 20/z~"^Mpc 
that are in the linear or mildly non-linear regime by the present time. In this paper we 
will only consider linear theory, in the form of the Zel'dovich approximation, to evolve 
the perturbations in the matter era. Matter accretion during the radiation era was not 
considered in the present paper. For linear perturbations the genus curve is isomorphic to 
the 'initial' genus curve which implies that the genus curve as a function of the number 
of standard deviations from mean density does not change with time if no additional 
defects enter the box. The effect of non-linear evolution on the genus curves is work in 
progress at the present moment (Avelino & Canavezes 1995). However, we expect this 
to be small on scales greater than 8h~^M.pc. We produce the density fluctuations in 
a comoving box whose size is chosen so that the genus curves obtained can be directly 
compared with those obtained from the CfA survey. We flx the scaling solution by flxing 
the length and number of defects of each type per horizon volume. These were put in the 
box with random positions and orientations. The number of defects per horizon volume 
is iV lb \/iV and the wavelength of the perturbations induced by them to be initially (at 
^eq) ^eq and |Aeq X |Aeq for the fllameuts and wakes respectively (Aeq = cteq). Structures 
seeded at later times will have a larger wavelength due to the scaling solution. These 
structures were produced by displacing the particles according to the symmetry point, 
axis and plane and with a dependence on the spatial coordinate given by each of the 
solutions of the Zel'dovich equation. In appendix 1 we give a more detailed description 
of how the fluctuations are produced for each toy model. In this article we assume the 
biasing parameter does not depend on the scale so that we can directly compare the genus 
curves obtained from the CfA survey with those obtained for the toy models considered 
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in this paper (at least for scales which are in the linear regime) which deal only with the 
distribution of the dark matter . 

4.4 Compensation 

Arbitrary energy momentum perturbations are not possible in a Robertson- Walker space- 
time (Trashen 1985; Trashen, Turok & Brandenberger 1986; Veeraghavan & Stebbins 
1990). When the strings are formed in the early universe the excess energy and linear 
momentum carried by the string is compensated by an equal deficit in the background 
radiation. The fact that energy and momentum must be conserved is well expressed by 
the following integral constraints 

/ 5p<N ^ 0, (42) 

I Sp^dN ?^ 0, (43) 
Jv 

where V is a volume much greater than the horizon volume. Because of the way in which 
we generate the perturbations, only positive density perturbations would be allowed if 
compensation was not taken into account. Although it is difficult to mimic the detailed 
effects of compensation in the genus curve we can at least, as a first approximation, cor- 
rect the genus curve for the average shift of the peak due to not including compensation 
by doing 

<^new = <^old - 4, (44) 

such that ((^new) = inside the box. 

5 Results 

5.1 Toy models and CDM 

In fig 3 we show a comparison of the density probability distribution for 3 of the models 
studied in the context of cold dark matter, averaged over 10 simulations. These were the 
5 filament model, the 5 wake model and the 25 sphere model. 

We can see from the pictures that for Ag = 6/z~"^Mpc the density distribution is very 
non- Gaussian and that it will approach Gaussianity as we go to larger smoothing scales. 
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We can also see that the largest departures from Gaussianity occur for the 5 filament 
model while for wakes the density probability distribution is very nearly Gaussian even 
at small scales. If all the models had the same number of objects per unit volume we 
would expect the sphere model to exhibit the largest departures from Gaussianity, while 
the wake model should be the most Gaussian. The reason for this is that while a point- 
like overdensity can only be inside or outside a chosen volume, line-like perturbations 
and planar perturbations can be partially inside several volumes (a planar perturbation 
with size 1x1 being able to touch more volumes than a line-like perturbation of size 
1). Consequently, the number of pieces contained in a chosen smoothing volume will be 
maximized for the wake model, making this model more Gaussian, and minimized for 
sphere model making it less Gaussian. The reason why our particular sphere model is 
more Gaussian than the filament model is that there are more objects per unit volume 
in the sphere toy model than in the filament toy model. 

The large departure from Gaussianity, especially on small scales, makes the genus 
curve very sensitive to the parametrization so that different genus curves are obtained if 
the genus is expressed directly as a function of the number of standard deviations from 
mean density or if we use the volume fraction parametrization for v. A comparison of 
two genus curves for the 5 filament model is shown in fig 4 in order to illustrate the effect 
of the parametrization by volume fraction. The genus curves are considerably different 
(the one parametrized by volume fraction being more like random phase), but as we have 
seen this effect should be most exaggerated for the 5 filament model at the smallest scale 
considered (Ag = 6/z~"^Mpc). For other toy models and smoothing scales the dependence 
of the genus curve on the parametrization will not be so considerable. However, for small 
smoothing lengths some information about the density probability distribution is lost. In 
this paper we shall use the volume fraction parametrization (unless it is said otherwise) 
used by Vogeley et al. (1994) when analysing the CfA observational data. 

In fig 5 and fig 6 we show a comparison of the genus curves, for our 6 CDM models of 
structure formation, with the predictions from the standard inflationary CDM model, as 
well as the observations taken from the CfA survey. Error-bars due to sample variance 
were not included for the sake of clarity. The sizes of the boxes for the toy models 
studied in this article were chosen such that a direct comparison between their genus 
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curves and CfA genus curves is possible for each of the smoothing scales considered. It is 
readily apparent from the graphs that from the models studied the one that best fits the 
observational data, at least from a topological point of view, is the 5 wake toy model. 
The amplitude of the genus curves seems to be only marginally larger than the CfA genus 
curves. Also, the width of the genus curves seems to mimic observations much better 
than the standard CDM model. The 10 wake model also does marginally better than 
standard CDM. The amplitude of the genus curves is always smaller than that predicted 
by the standard CDM scenario and in better agreement with observations. 

The standard CDM model gives too large amplitude of the genus curve, especially 
on small scales and it fails to match other features of the observed genus curves. For 
example, the standard CDM model gives Wi^ ~ 2 and although sample variance allows 
some fluctuations around this value it is not enough to explain genus peak widths as 
large as 2.5 or 2.6 with a very large confldence level (> 90%) (Vogeley et al. 1994). It 
is possible to have other random phase models with a smaller amplitude of the genus 
curve, as in open models or models with a non-zero cosmological constant. What hap- 
pens in this case is that the growth of density perturbations in the linear regime slows 
down while the non-linear growth of perturbations proceeds in the normal way. So, for 
the same normalization of the spectrum of density perturbations non-linear perturba- 
tions are more non-linear in these models. Non-linearities introduce correlations between 
phases of different Fourier modes, decreasing the number of independent structures and 
consequently reducing the genus amplitude. However, as found by Vogeley et al. (1994) 
the problem of matching the other statistics, especially the width of the genus peak, 
remains unsolved. 

The 5 and 10 fllament models do better than spheres and marginally worse than 
standard CDM. The genus amplitude is higher than for standard CDM. However, these 
models provide a better fltting to the width of the observed genus curves than the stan- 
dard CDM inflationary scenario. The sphere model is clearly ruled out. The amplitude 
of the genus curves is too large and it fails to match the shape of the observational genus 
curves. 

Also apparent from the graphs is that for most of the smoothing scales considered, 
the genus amplitude is an increasing function of the number of defects. This should be 
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expected because we are increasing the number of structures present inside the box. 

In fig 7 we plotted the genus curves obtained for the best model (the 5 wake model) 
now with the error bars properly included. The line represents the average genus curve 
among 10 realizations of the model. The error bars are one sigma error bars over these 
realizations. The asymmetry of the genus curves is visible, as we would expect from 
almost any non-Gaussian model of structure formation. However, the parametrization 
by volume fraction makes the genus curves look less asymmetric, especially on small 
scales, making the curves more like random-phase models. 

In fig 8 we show a statistical comparison of the genus curves for the wake model with 
the predictions from the CfA survey. The genus amplitude for this model is marginally 
larger than the observed genus amplitude for some the scales considered, although it does 
much better than standard CDM at matching the amplitude of the CfA genus curves. 
Again, the error bars we see in the plots are one sigma error bars of 10 simulations. 
The shift of the genus curves does not seem to be in agreement with observations for 
most of the scales considered. The width of the genus peak seems to be consistent with 
observations from the CfA survey for all of the scales considered but in this case the 
error bars due to sample variance are very large so that a very wide range of genus peak 
widths are possible. This toy model fits the observations better than any random phase 
model tried by Vogeley et al. (1994). 

5.2 Toymodels and HDM 

In fig 9 we show a comparison of the density probability distribution for 3 of the models 
studied as a function of the smoothing length and averaged over 10 simulations. These 
were the 5 filament model, the 5 wake model and the 25 sphere model with hot dark 
matter. These models are more Gaussian than the corresponding models with CDM and 
consequently the dependence of the genus curves on the parametrization is smaller in 
this case. For the wake model the density probability distribution is nearly Gaussian 
even for Ae = 6h~^M.pc. 

In fig 10 and fig 11 we show a comparison of the genus curves for 6 models of structure 
formation we studied in the context of hot dark matter with the predictions from the 
standard CDM model and observations taken from the CfA survey. The amplitude of the 
genus curves is smaller with HDM than with CDM, especially on small scales, because 
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an adittional smoothing was introduced due to the free-streaming of the neutrinos. The 
models that best mimic the observational results are again the wake models (particularly 
the 10 wake model). 

In fig 12 we can see the genus curves obtained for the best model (the 10 wake model) 
now with the error bars properly included. These seem to be more like random-phase 
than the genus curves obtained for the 5 wake model. The line represents the average 
genus curve among 10 realizations of the model and the error bars are one sigma error 
bars. 

Figure 11 shows a statistical comparison of the genus curves for the 10 wake model 
with the predictions from the CfA survey. The genus amplitude for this model seems to 
be consistent with the observed genus amplitude for most of the scales considered. Again, 
the error bars we see in the plots are one sigma error bars of 10 simulations. The shift 
of the genus curves does not seem to be in agreement with observations for most of the 
scales considered. The width of the genus peak seems to be consistent with observations 
from the CfA survey although in this case the error bars due to sample variance are very 
large. This toy model also fits the observations better than any random phase model 
tried by Vogeley et al. (1994). 

6 Discussion 

We can conclude from the results presented in this article that the genus statistic is a 
good discriminator between Gaussian and non-Gaussian models of structure formation. 
In addition to this, it was shown to be a good statistic to distinguish between different 
toy models of structure formation, sensitive to the shape, number of structures seeded, 
and dark matter type. We also showed that at least for some of the models considered (in 
particular the wake models) that there is a better agreement with the observations than 
that verified for random phase models of structure formation. It is necessary to use string 
network simulations in combination with numerical codes that generate and evolve the 
density perturbations seeded by such networks in order to properly test the cosmic string 
model for structure formation. However, we have shown that there are some features on 
the observed genus curves (e.g. genus peak width considerably larger than 2 and genus 
amplitude smaller than standard CDM), that cannot be easily matched by random phase 
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models of structure formation, but which are matched by some toy models considered 
in this paper over a range of smoothing lengths. Although this cannot be considered 
to be a serious quantitative test of the cosmic string paradigm for structure formation, 
it provides a good motivation for arguing that some statistical features of topological 
defect models of structure formation are better than random phase models predicted by 
most inflationary scenarios. 
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Appendix 1: Generation of perturbations for the toy models 



In the sphere model a number N ± y N of spheres were put, per horizon per Hubble 
time, at random positions in the box. For each sphere we chose a position p at random. 
Particles were then displaced according to the Zel'dovich approximation. Consider a 
sphere laid down at an instant t. A particle at a position p' would move from that time 
t till the present time to a comoving distance given (in linear theory) by 

^ = ^GMt,2^6,(t), (45) 
|x| 

where x = p' — p and 

^ 2 t 3 fto^^'^' 



If we account for the growth of the mass of the loops chopped off by the network of 
cosmic strings (M oc t) we have that the quotient of the perturbations laid down at two 
different times ti and ^2 at the same position in space is given by 

^(t2) hxhs{t2)' 

We can see that in the case of the sphere model (only) perturbations seeded at later 
times can have larger amplitude. 

In the fllament model for each fllament a point p and a unit vector v were chosen at 
random. Particles were then displaced according to the Zel'dovich approximation. Let 
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us consider a particle at a position p' and let us define the vector y = p' — P- Consider 
the vector x = y — (y • v)v and assume that the filament has a size given by Sf. If 
|(y ■ v)v| < 0.55"/ then a particle at a position p' would move a comoving distance given 
by 

6 , .2 X 



ij = --GMLU'—^bf{t). (48) 



5 



X 



where 



\ 2/3 \ Q / J. \ 2/3 

to \ 1 MO \ o I i^O 




If |(y • v)v| > 0.55"/ the particle would not move at all. The mass per unit length of the 
cosmic strings is approximately constant over time {Ml = const) and so we have that 
the quotient of the perturbations laid down at two different times ti and t2 at the same 
position in space is given by 

= MM. (50) 

Consequently, the perturbations seeded at earlier times will have larger amplitude. The 
initial filament comoving size (at teq) was taken to be Sf = Aeq (where Aeq = cteq) and 
its size increases with time proportionally to the horizon so that Sf oc t^/'^ . So, at later 
times larger structures are formed but those will be less dense. 

In the wake model for each filament a point p and a unit vector v were also chosen 
at random. Particles were then displaced according to the Zel'dovich approximation. 





Let us consider the vector A = ( | and a vector perpendicular to it B 
Consider the rotation matrices 



M2 = I cos^ sin^ I , (51) 
— sin^ cos ^ 

cos l3 sin [i 
Ml = I 1 I . (52) 
■ sin l3 cos [i 

The vector A' = M2 • Ml • A obtained as the result of multiplying the matrices Ml 
and M2 by A is a new vector given by A' = (sin /3, sin ^ cos /3, cos ^ cos /3). The vector 
B' = M2 • Ml • B is perpendicular to A' and is given by 

cos l3 sin a 

B' = I cos ^ cos a — sin ^ sin /3 sin a | . (53) 
sin 6 cos a — cos 6 sin [i sin a 
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The angles (3 and 6 were chosen subject to the constraint v = A' and a was choosen at 
random in the interval < a < 2tt. We have now two perpendicular vectors A' and B' 
and we need to find C perpendicular to this two such that 

C'.A' = 0, (54) 



C'.B' = 0, (55) 

and 

|C'| = 1. (56) 

To see how a particle at a position p' will move let us consider the vector y = p' — p 
and the vectors Xi = y • B', X2 = y • C and x = y — (y • Xi)xi — (y • X2)x2. Let us 
assume the size of the wake to be Sw If |xi| < 0.53^ and |x2| < 0.53^ then a particle 
at a position p' would move a comoving distance given by 

^ 2 X 

tp = --u^t^ — hyj{t). (57) 
5 |x| 

where 

^■"")=((^)"'-l)- <^**) 

If |xi| > Q.^3w or |x2| > 0.5S'„, the particle would not move at all. The quotient of the 
perturbations laid down at two different times ti and ^2 at the same position in space is 
given by 

We can see again that the perturbations seeded at earlier times will have larger amplitude. 
The initial wake size was taken to be |Aeq and it grows proportionally to the horizon as 
in the filament case. 

Appendix 2: The simulation boxes 

As we have said before the size of the simulation boxes was chosen in a way to enable 
direct comparison with the results from Vogeley et al. (1994) for the CfA survey. To 
ensure that the topology was not dominated by shot noise the smoothing length must 
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Table 1: Volume statistics for the CfA survey (note that the total volume and consequently the 
number of resolution elements for the toy model simulations is the same as for the CfA survey). 

be larger than the average intergalaxy (or interparticle) separation. To determine the 
maximum distance appropriate for a given choice of smoothing length Vogeley et al. 
(1994) found 

'"max such that 

?^(?"max) = Ae~^. (60) 

This means that the number of galaxies included is an increasing function of the smooth- 
ing length. Table 1 shows the volume of the survey as a function of the smoothing length. 
The number of resolution elements, defined by iVj-es = ~%Tr^, and the number of galax- 
ies included in the topological analysis of the CfA survey are also given as a function 
of the smoothing length. The volume of the simulation boxes and so the number of 
resolution elements were chosen to be the same as for the CfA survey analysis so that 
direct comparison between our results and observations was possible. The smoothing 
length was always more than 2 times larger than the average interparticle spacing. 
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Fig. 1 Comparison of the genus curve obtained for a 64^ simulation of a P{k) cx k 
power spectrum smoothed on a scale Ag equal to \/2 times the grid spacing, obtained 
(a) analytically (b) numerically using our program to calculate the genus. 
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Fig. 2 Isodensity contours with S = 1.5i^ for (a) the 5 filament model b) the 5 
wake model c) the 25 sphere model with cold dark matter. The box size is 3.13 x 
10®/i~^Mpc^, and the smoothing length is Ag = 6/i~^Mpc. 
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Fig. 3 Density probability distribution for several models studied with cold dark 
matter as a function of v calculated directly from the variance of the density distri- 
bution (a) 5 filament model (b) 5 wake model and (c) 25 sphere model. 
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Fig. 4 Comparison of genus curves for the 5 filament model assuming the dark 
matter to be cold, (a) using the volume fraction to prametrize v (b) Calculating v 
directly from the variance of the density distribution. 
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Fig. 5 Comparison of genus curves for the 5 wake, 5 filament and 25 sphere models 
with both standard CDM genus curves and genus curves obtained from the CfA 
survey. Error bars are not included for clarity (see fig 8). It is clear from the picture 
that the most favoured model is the 5 wake model. 
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Fig. 6 Comparison of genus curves for the 10 wakes, 10 filaments and 50 spheres 
models with both standard CDM genus curves and genus curves obtained from the 
CfA survey. 
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Fig. 7 Genus curves for the 5 wake model. Line represents the average genus curve 
among 10 realizations of the model. Error bars are one-sigma. 
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Fig. 8 Statistical comparison of genus curves obtained from the CfA survey with 
genus curves obtained from 10 realizations of the 5 wake model. Error bars on the 5 
wakes model are one-sigma. 
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Fig. 9 Density probability distribution for several models studied with hot dark 
matter as a function of v calculated directly from the density distribution (a) 5 
filament model (b) 5 wake model and (c) 25 sphere model. 
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Fig. 10 Comparison of genus curves for the 5 wake, 5 filament and 25 sphere 
models in the context of hot dark matter with both standard CDM genus curves and 
genus curves obtained from the CfA survey. Error bars are not included for clarity 
(see fig 13). 
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Fig. 11 Comparison of genus curves for the 10 wake, 10 filament and 50 sphere 
models in the context of hot dark matter with both standard CDM genus curves and 
genus curves obtained from the CfA survey. It is clear from the picture that the most 
favoured model is the 10 wake model. 

34 




Fig. 12 Genus curves for the 10 wake model with HDM. Line represents the average 
genus curve among 10 realizations of the model. Error bars are one-sigma. 
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Fig. 13 Statistical comparison of genus curves obtained from the CfA survey with 
genus curves obtained from 10 realizations of the 10 wake model with HDM. Error 
bars are one-sigma 
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